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Tikuisis  P,  Nishi  RY,  Weathersby  PK.  Use  of  the  maximum  likelihood  method  in  the  analysis 
of  chamber  air  dives.  Undersea  Biomed  Res  1988;  15(4):301-313.^The  method  of  maximum 
likelihood  was  used  to  evaluate  the  risk  of  decompression  sickness  (DCS)  for  selected  chamber 
air  dives.  The  parameters  of  two  mathematical  models  for  predicting  DCS  were  optimized 
until  the  best  agreement  (as  measured  by  maximum  likelihood)  corresponding  to  the  observed 
DCS  incidents  from  a  series  of  dives  was  attained.  The  decompression  data  used  consisted  of 
800  man-dives  with  21  incidents  of  DCS  and  6  occurrences  of  marginal  symptoms.  The  first 
model  investigated  was  based  on  a  nonlinear  gas  exchange  in  a  series  arrangement  of  four 
compartments.  The  second  model  was  based  on  a  monoexponential  gas  exchange  in  a  parallel 
arrangement  of  two  compartments.  The  overall  statistical  success  in  describing  the  800  man- 
dives  was  quite  similar  for  the  two  models.  Predictions  of  safety  for  dives  not  part  of  the 
original  data  differed  for  the  models  due  to  differences  in  gas  kinetics.  For  short,  no¬ 
decompression  dives,  the  series  arrangement  of  compartments  predicted  a  lower  incidence  of 
DCS.  These  predictions  were  more  consistent  with  the  outcome  of  subsequent  testing  than 
were  predictions  of  the  parallel  compartment  model.  Predictions  of  the  series  arrangement 
model  were  also  similar  to  those  of  a  single-compartment,  two-exponential  model  that  was 
evaluated  with  over  1700  man-dives  by  the  U.S.  Navy. 


decompression  sickness 
dive 

maximum  likelihood 


models 

tables 

probability  Q,  f . 


The  variability  in  the  occurrence  of  decompression  sickness  (DCS)  makes  the 
prediction  of  DCS  difficult  for  any  given  dive  profile.  Conventional  methods  for 
establishing  the  risk  of  DCS  for  a  particular  dive  require  a  great  number  of  man- 
dives.  For  example,  if  in  100  man-dives  for  a  particular  trial,  5  incidents  of  DCS 
occurred  (5%  occurrence),  then  the  95%  confidence  interval  limits  for  this  occurrence 
indicate  that  the  probability  of  DCS  can  be  between  1 .64  and  1 1 .28%  ( I ).  To  “tighten” 
the  prediction  of  the  incidence  of  DCS  for  this  particular  dive  trial,  many  more  man- 
dives  would  be  required.  For  example,  65  incidents  of  DCS  out  of  1300  man-dives 
(5%  occurrence)  would  narrow  the  95%  confidence  interval  to  between  4  and  6%  for 
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a  particular  dive  trial.  Since  most  dive  trials  involve  considerably  fewer  than  100 
individuals,  a  method  for  combining  dive  trials  in  a  single  analysis  is  required. 

Weathersby  et  al.  (2)  suggested  using  the  method  of  maximum  likelihood  (3)  to 
evaluate  the  risk  of  DCS.  This  method  applies  a  predictive  model  for  DCS  to  actual 
outcomes  from  a  series  of  dive  trials.  Inputs  to  the  predictive  model  are  a  decompres¬ 
sion  dose  response  function  and  any  decompression  model  that  can  be  stated  math¬ 
ematically.  The  decompression  dose  response  function  is  a  measure  of  the  cumulative 
risk  of  DCS  for  a  dive.  The  parameters  of  the  predictive  model  are  adjusted  until  the 
best  agreement  or  maximum  likelihood  corresponding  to  the  DCS  events  of  past  trials 
is  obtained.  Once  the  parameters  of  the  predictive  model  have  been  optimized  using 
any  number  of  dives  from  different  sources,  the  model  can  be  applied  to  predict  the 
probability  of  DCS  for  any  new  (tested  and  untested)  dive  profile. 

In  addition,  different  decompression  models  can  be  tested  to  determine  which 
model  best  fits  the  data.  In  this  study,  the  method  of  maximum  likelihood  is  used  to 
compare  the  prediction  of  DCS  between  a  nonlinear  series  compartments  (NLSC) 
model  (4,  5),  featuring  four  compartments  in  series,  and  a  monoexponential  parallel 
compartments  (MEPC)  model,  used  by  Weathersby  et  al.  (6),  featuring  two  com¬ 
partments  in  parallel.  Although  many  other  models  involving  more  than  two  com¬ 
partments  in  parallel  exist  (4),  they  offer  little  or  no  statistical  advantage  over  the 
two-compartment  MEPC  model  in  the  present  analysis  (6),  but  may  be  justified  by 
larger  data  sets  (7).  The  method  of  maximum  likelihood  has  not  been  previously 
applied  to  compartments  in  series  such  as  the  NLSC  model.  The  decompression 
model  that  gives  the  largest  maximum  likelihood  value  is  defined  as  the  one  that  best 
fits  the  data.  For  specific  classes  of  models,  formal  statistical  tests  allow  rigorous 
comparison  of  model  predictions  (2, 3).  A  desirable  feature  of  the  method  of  maximum 
likelihood  is  that  any  set  of  dive  data  can  be  used  for  the  comparison  among  models. 
For  instance,  the  MEPC  model  can  be  tested  against  dive  trials  that  followed  the 
decompression  procedures  governed  by  the  gas  kinetics  of  the  NLSC  model,  and 
vice  versa. 

In  summary,  a  probabilistic  model  optimized  by  maximum  likelihood  can  be  used 
to  test  decompression  procedures  and  to  predict  the  probability  of  DCS.  It  can  be 
used  as  an  optimizing  tool  for  determining  the  parameters  of  decompression  models 
and  as  an  evaluative  tool,  but  the  prediction  must  be  tested  and  the  actual  incidence 
must  be  established  by  experiment.  This  study  examines  the  use  of  the  method  of 
maximum  likelihood  as  applied  to  the  NLSC  and  MEPC  models  in  optimizing  model 
parameters,  evaluating  dive  profiles,  and  comparing  the  prediction  of  DCS  for  par¬ 
ticular  dive  tables. 


LIKELIHOOD  CRITERION 

The  variability  in  DCS  can  be  described  in  terms  of  probabilities.  The  probability 
of  DCS  is  a  function  of  the  dive  profile  and  the  decompression  model  tested,  i.e.. 
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1  —  P  DCS-  (2) 

If  Y  is  the  actual  outcome  of  DCS  for  a  dive,  then  the  predicted  probability  of  the 
actual  outcome  or  observation  is 


P outcome  —  (P Dcs)^  (1  —  P DCs)'  ' 


(3) 


where  Y  =  0  if  DCS  does  not  occur,  and  Y  =  I  if  DCS  does  occur.  For  those  dives 
that  defy  a  definitive  diagnosis  (sometimes  called  marginal  symptoms  or  niggles)  an 
assignment  of  Y  =  0.5  can  be  made.  In  other  analyses  such  a  choice  has  not  been 
shown  to  affect  the  ultimate  conclusions  of  the  analysis  (2,  8). 

Assuming  that  the  outcome  of  each  dive  is  independent,  the  total  predicted  prob¬ 
ability  of  all  the  known  outcomes  is  given  by  the  likelihood  function, 


l  =  nPo, 


(4) 


where  n  is  the  total  number  of  man-dives.  L  is  usually  a  very  small  number,  but 
greater  than  zero,  because  it  is  the  product  of  many  numbers  less  than  1.  It  is 
customary  to  use  the  natural  logarithm  of  L,  or  LL  (log  likelihood). 


If  ^ In  P outcome*  (3) 

Since  L  is  less  than  1,  LL  is  negative.  The  maximum  likelihood  value,  which  defines 
the  best  fit  between  the  model  and  data,  is  given  by  the  smallest  absolute  value  of  LL 
(or  smallest  negative  value)  as  shown  in  Fig.  1.  The  modified  Marquardt  algorithm 
(9)  was  used  to  obtain  this  value  by  varying  the  parameters  of  the  predictive  model 
until  the  best  agreement  with  the  data  was  found. 


Risk  Coefficient,  104C 

Fig.  I.  Log  maximum  likelihood.  LL.  plotted  against  the  risk  coefficient.  C.  for  the  NLSC  model. 
Solid  circles  indicate  the  95%  confidence  limits  based  on  the  Likelihood  Ratio  Test  0). 
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To  use  the  method  of  maximum  likelihood,  the  decompression  dose  response 
function  or  PXs  (see  Eq.  /)  must  be  defined.  The  function  that  was  adopted  for  this 
study  is  known  as  the  “risk”  function  (2,  6,  10), 


where  r  is  the  measure  of  risk  and  t  is  time.  If  r  is  large,  then  over  a  period  of  time 
the  probability  of  DCS  approaches  1.  For  each  decompression  model  examined,  r  is 
defined  empirically  and  its  value  is  dependent  on  the  actual  decompression  profile 
(see  Decompression  Models).  Generally,  r  is  non-zero  when  a  state  of  inert  gas 
supersaturation  occurs.  Depending  on  the  degree  and  length  of  this  state,  PD cs  will 
increase  accordingly.  The  calculation  of  /Vs  is  continued  until  r  falls  below  zero,  at 
which  point  the  value  of  Foes  represents  the  predicted  cumulative  probability  of  DCS 
for  the  dive. 

For  models  that  can  be  expressed  as  simplified  or  more  general  versions  of  other 
models,  a  formal  Likelihood  Ratio  Test  can  establish  whether  a  model  is  a  significant 
improvement  over  another  in  describing  a  data  set  (2,  3).  For  example,  to  justify 
increasing  the  generality  of  a  model  by  estimating  one,  two,  three,  etc.,  additional 
parameters,  the  LL  must  increase  by  at  least  1 .92,  3.00, 3.91  units,  etc.,  respectively. 


DECOMPRESSION  MODELS 


1.  Nonlinear  series  compartments  model 

The  NLSC  model  for  air  diving  (4,  5)  consists  of  four  compartments  in  series 
representing  tissues  having  different  rates  of  gas  exchange.  The  rate  of  change  of  gas 
pressure  in  the  four  compartments  is  described  by  the  following  set  of  4  nonlinear 


differential  equations: 

dPJdt  =  A[(B  +  P0  +  P i)(F»  -  Pi)  -  (B  +  Pi  +  P2)(Pi  -  P2)],  (7) 

dPJdt  =  A[(B  +  P,  +  P2)(P,  -  P2)  -  (B  +  P2  +  P,)(P2  -  P,)l  (8) 

dPJdt  =  A((B  +  P2  +  Py)(P2  -  Py)  -  (B  +  P,  +  P«)(P,  -  P«)],  (9) 

dPJdt  =  A[(B  +  P,  +  P«)(P3  -  P«),  (10) 


where  P(  is  the  total  gas  pressure  in  compartment  i,  P0  is  the  current  ambient  pressure, 
and  A  and  B  are  constants  (0.0026  min-1  •  ATA-1  and  8.31  ATA,  respectively). 
Equations  7-10  are  those  that  control  the  uptake  and  elimination  of  gas  for  the  original 
Kidd-Stubbs  pneumatic  analogue  decompression  computer  model  (4),  the  Kidd- 
Stubbs  1971  decompression  tables  (11),  and  the  current  Defence  and  Civil  Institute 
of  Environmental  Medicine  (DCIEM)  1983  decompression  tables  (12).  The  only 
differences  among  these  3  cases  are  in  the  safe  ascent  criterion  for  each  compartment. 
The  measure  of  risk  for  the  present  likelihood  analysis  is  defined  as 

r  =  r,  +  r2  +  r,  +  r„  (11) 

where,  following  the  method  of  Weathersby  et  al.  (6),  the  risk  for  compartment  /  is 
assumed  to  be 
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r,  =  C(0. 79P,  -  P0)/P0,  (12) 

and  where  C  is  the  risk  coefficient  to  be  estimated.  The  risk  defined  by  Eq.  12  is 
proportional  to  the  relative  supersaturation  of  the  N2  component,  hence  the  factor 
0.79.  If  r,  <  0,  then  r,  is  set  to  0.  The  value  of  C  is  estimated  by  fitting  the  model  to 
the  actual  time-depth  profile  data  for  each  dive  until  the  maximum  value  of  LL  is 
obtained.  The  model  constants,  A  and  B,  can  also  be  optimized  using  the  maximum 
likelihood  method;  however,  time  constraints  did  not  allow  this. 

Figure  2  graphically  demonstrates  the  NLSC  model  with  C  =  0.000154  •  min-1  to 
a  single  dive.  This  dive,  to  185  fsw  (5.60  ATA)  with  a  bottom  time  of  125  min,  had  4 
divers  and  no  incidence  of  DCS.  At  the  end  of  decompression  (/  =  827  min),  the 
predicted  incidence  of  DCS  is  5.97%,  while  r  is  still  greater  than  zero  owing  to  the 
supersaturation  of  N2  in  the  model  compartments.  After  1 100  min,  r  according  to  Eq. 
12,  becomes  zero  and  the  final  predicted  incidence  of  DCS  is  7.04  ±  1.55%  (95% 
confidence  limit). 

2.  Mono-exponential,  parallel  compartments  model 

The  MEPC  model  is  based  on  an  arrangement  of  single  exponential  compartments 
in  parallel.  Weathersby  et  al.  (6)  found  that  two  compartments  were  adequate  to  fit 
the  present  data  in  their  study.  Hence,  in  this  study  a  two-compartment  mono- 
exponential  model  [model  3  in  (6)]  is  used.  The  exchange  of  inert  gas  in  each  of  the 
two  compartments  is  described  by  the  following  linear  differential  equation: 


Time  (min) 


Fig.  2.  Measure  of  risk  of  DCS.  r.  and  the  probability  of  incidence  of  DCS.  Foes,  predicted  by  the 
NLSC  model  for  a  dive  to  I8J  fsw  for  125  min  plotted  as  a  function  of  time.  The  time-depth  profile  of  the 
dive  and  the  average  gas  pressure  of  the  four  compartments  of  the  NLSC  model  are  also  shown. 
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dPsJdt  =  -(/>*,  -  0J9Po)lk„  (13) 

where  PN,  is  the  pressure  of  the  N2  component  of  gas  in  compartment  i  and  kt  is  the 
time  constant  of  the  N,  uptake  and  elimination  of  compartment  /. 

The  measure  of  risk  is  defined  as 


r  =  r,  +  r2  (14) 

where 

r,  =  C,(7V  -  Wo,  (15) 

and  C j  is  the  risk  coefficient  of  compartment  i.  The  values  of  C ,  and  k,  were  estimated 
by  fitting  the  model  to  the  data  until  the  maximum  value  of  LL  was  obtained. 

RESULTS 

The  dive  data  were  obtained  from  the  CANDID  (Canadian  Dive  Data)  data  bank 
(13)  and  are  the  same  data  used  in  data  set  “D”  in  (6).  The  data  used  consist  of  163 
air  dive  profiles  in  a  hyperbaric  chamber  totaling  800  man-dives  that  were  conducted 
during  1967-1968  with  the  Kidd-Stubbs  pneumatic  analogue  decompression  computer 
and  other  decompression  criteria  (4).  The  dive  depths  for  this  data  set  ranged  from 
99  to  300  fsw  (3.00  to  9.08  ATA)  and  the  bottom  times  ranged  from  6  to  360  min. 
There  were  21  incidents  of  DCS  and  6  occurrences  of  marginal  symptoms  of  DCS. 
Incidents  of  DCS  were  assigned  an  outcome  of  Y  =  1.0  and  marginal  cases  were 
assigned  an  outcome  of  Y  -  0.5. 

Using  these  data,  the  values  of  LL  are  plotted  against  the  risk  coefficient.  C,  for 
the  NLSC  model  in  Fig.  1.  The  maximum  value  of  LL  (=  - 102.57)  was  obtained 
with  C  =  0.000154  -  min'1.  This  value  was  used  to  generate  Fig.  2.  For  the  MEPC 
mode],  the  maximum  value  of  LL  ( =  -  100.63)  was  obtained  with  k,  =  3.91  min,  k2 
=  382  min,  C,  =  0.00615  •  min'1,  and  C2  =  0.00126  ■  min'1  (6).  The  values  of 
likelihood  for  the  2  models  are  close  enough  that  the  additional  three  estimated 
parameters  of  the  MEPC  model  would  not  be  justified  by  the  Likelihood  Ratio  Test 
if  it  could  be  applied.  Since  the  models  are  not  suitable  for  this  test,  then  within 
current  limits  of  statistical  theory,  the  results  can  only  be  declared  as  similar. 

It  is  interesting  to  compare  the  predictions  of  these  2  models  under  different 
categories  of  risk  of  DCS.  Figure  3  shows  the  results  of  the  2  models  after  separating 
the  predictions  of  DCS  (using  Eq.  6)  into  four  categories:  <2%,  2-5%,  5-10%,  and 
>10%.  For  example,  the  NLSC  model  predicted  a  DCS  incidence  of  less  than  2% 
for  312  out  of  the  800  man-dives  in  the  data  set.  Of  these  312  dives,  the  average 
predicted  incidence  was  1.47%  and  the  average  observed  incidence  was  1.44%.  In 
comparison,  the  MEPC  model  predicted  a  DCS  incidence  of  less  than  2%  for  238  out 
of  the  800  man-dives,  and  of  these  the  average  predicted  incidence  was  1 .36%  while 
the  average  observed  incidence  was  0.84%.  No  dives  were  predicted  to  have  a  DCS 
incidence  of  >10%  using  the  NLSC  model  in  contrast  to  the  15  man-dives  predicted 
to  exceed  10%  DCS  incidence  by  the  MEPC  model.  Despite  these  differences,  the 
chi-square  goodness-of-fit  test  for  agreement  between  model  prediction  and  actual 
outcome  for  all  4  categories  is  0.60  for  the  NLSC  model  and  1.29  for  the  MEPC 
model;  these  low  values  indicate  that  both  models  predict  outcome  of  DCS  equally 
well. 
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Fig.  3.  Average  observed  incidence  of  DCS  plotted  against  the  average  predicted  incidence  of  DCS 
for  the  NLSC  model  (open  circles)  and  MEPC  model  (solid  circles)  in  four  categories:  <2%,  2-5%,  5- 
10%,  and  >10%.  Broken  line  shows  the  line  of  agreement  and  vertical  lines  indicate  the  95%  confidence 
interval  on  the  observed  binary  outcome. 


Application  of  Theory 

Once  the  parameters  of  the  predictive  model  have  been  optimized,  the  model  can 
be  used  to  predict  the  probability  of  DCS  for  a  series  of  different  but  related  dives. 
Figure  4  shows  the  probability  of  DCS  predicted  by  both  the  NLSC  and  MEPC 
models  for  the  DCIEM  no-decompression  air  diving  limits  (12)  as  shown  in  Table  I. 
These  no-decompression  limits  are  more  conservative  than  those  from  the  USN 
Standard  Air  Decompression  Tables  (14,  15)  (also  shown  in  Table  1).  For  depths  of 
30  to  40  fsw  (0.91  and  1.21  AT A),  which  have  not  been  tested,  the  predicted  probability 
of  DCS  is  noticeably  higher  than  that  for  deeper  depths;  hence,  according  to  this 
analysis  the  no-decompression  limits  (12)  from  these  shallow  depths  entail  a  greater 
risk  of  DCS. 

Dives  deeper  than  50  fsw  (1.51  ATA)  have  been  tested  according  to  the  DCIEM 
no-decompression  air  diving  limits.  No  instances  of  DCS  occurred  in  the  131  man- 
dives  considered  here  (16),  an  outcome  that  has  95%  confidence  limits  of  0  to  2.8% 
DCS  (1).  In  this  depth  range  the  NLSC  model  predicts  a  uniformly  low  probability 
of  DCS  averaging  0.2%,  which  is  quite  consistent  with  the  outcome.  This  behavior 
is  not  too  surprising  since  the  dive  limits  were  obtained  by  a  nonprobabilistic  model 
with  the  same  kinetic  constants  as  the  present  NLSC  model  (5,  12).  On  the  other 
hand,  the  MEPC  model  predicts  a  logarithmically  increasing  probability  of  DCS  with 
depth  that  averages  3.1%,  or  slightly  outside  the  observed  outcome.  The  reason  for 
the  contrasting  trend  between  the  2  models  lies  with  the  difference  in  the  kinetics  of 
N2  gas  uptake  and  elimination  in  the  compartments  of  the  2  models.  The  series 
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TABLE  1 

Comparison  of  No-Decompression  Air  Diving  Limits 


Depth,  fsw 

Maximum  Bottom  Time,  min 

DCIEM* 

USN 

30 

380 

_ 

40 

175 

200 

50 

75 

100 

60 

50 

60 

70 

35 

50 

80 

25 

40 

90 

20 

30 

100 

15 

25 

110 

12 

20 

120 

10 

15 

130 

8 

10 

140 

7 

10 

150 

7 

5 

160 

6 

5 

170 

5 

5 

180 

5 

5 

•These  dive  limits  have  been  tested  from  50  to  180  fsw  at  DCIEM. 


compartmental  structure  of  the  NLSC  model  leads  to  a  slower  rate  of  gas  uptake  in 
the  compartments,  whereas  the  “fast”  compartment  ( k ,  =  3.91  min)  of  the  MEPC 
model  saturates  quickly.  It  is  this  fast  compartment  that  increases  ‘he  calculated 
measure  of  risk  substantially  more  than  do  the  compartments  of  the  NLSC  model. 

Figure  4  also  shows  the  predictions  of  model  5  in  (6),  a  single-compartment,  two- 
exponential  model  evaluated  with  over  1700  man-dives  used  to  estimate  the  risk  of 
bends  for  different  air  decompression  tables  (15).  In  the  depth  range  of  greater  than 
50  fsw  (1.51  ATA),  the  average  prediction  is  1.2%  DCS  which,  along  with  the  NLSC 
model,  is  consistent  with  the  outcome  for  the  present  test  dives. 

Differences  in  model  predictions  are  not  as  pronounced  for  dives  with  decompres¬ 
sion  stops  because  the  slower  compartments  predominate.  Dives  with  decompression 
stops  are  usually  long  enough  for  the  gas  pressure  in  the  compartments  of  the  NLSC 
model  to  increase  the  measure  of  risk  to  a  degree  similar  to  that  predicted  by  the 
MEPC  model.  For  example,  of  the  800  man-dives  examined  earlier,  the  average 
predicted  incidence  of  DCS  for  the  NLSC  model  was  2.92%  and  that  for  the  MEPC 
model  was  2.97%.  However,  there  are  still  some  important  differences  between  the 
predictions  of  the  2  models  for  dives  with  decompression  stops. 

This  is  demonstrated  in  Fig.  5,  where  the  probability  of  DCS  predicted  by  both  the 
NLSC  and  MEPC  models  for  standard  air  dives  to  150  fsw  (4.54  ATA)  from  the 
DCIEM  1983  decompression  tables  (12)  is  shown  as  a  function  of  bottom  time  at  150 
fsw.  The  NLSC  model  predicts  a  monotonically  increasing  probability  of  DCS  with 
increasing  bottom  time.  The  MEPC  model,  on  the  other  hand,  predicts  a  high  prob¬ 
ability  of  DCS  for  a  short  bottom  time  of  7  min,  then  decreasing  probabilities  until  a 
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Bottom  Time  (min) 

Fig.  4.  Probability  of  the  incidence  of  DCS  plotted  against  depth  for  the  DCIEM  no-decompression 
air  diving  limits  showing  the  difference  between  the  NLSC  model  ( solid  line),  the  MEPC  model  ( broken 
line),  and  model  5  in  (6)  (dolled line)  predictions. 

minimum  is  reached  at  a  bottom  time  of  25  min,  followed  by  increasing  probabilities 
thereafter.  In  all  cases,  the  MEPC  model  predictions  of  DCS  are  higher  than  those 
for  the  NLSC  model.  It  is  interesting  to  note  that  the  two  model  predictions  are  in 
close  agreement  in  the  mid-range  of  bottom  time,  say  25  to  55  min,  which  is  the  region 
where  data  were  available  (in  the  data,  dives  between  145  and  155  fsw  averaged  30 
min  in  bottom  time)  but  diverge  outside  this  range.  This  difference  results  from  the 
different  gas  kinetics  of  the  2  models  discussed  above.  Actual  testing  of  these  profiles 
(17)  produced  1  case  of  DCS  out  of  20  man-dives  with  20  and  30  min  bottom  times 
(95%  confidence  limits  of  0. 1  to  25%  DCS).  Both  NLSC  and  MEPC  model  predictions 
are  consistent  with  that  outcome. 

Also  shown  in  Fig.  5  are  the  predictions  (15)  for  model  5  in  (6),  which  are  also 
consistent  with  the  outcome  cited  above  (17).  These  values  follow  a  similar  trend  to 
the  MEPC  model  prediction  but  with  much  lower  values  around  20  min  bottom  time 
and  progressively  much  higher  values  for  bottom  times  beyond  35  min.  These  differ¬ 
ences  are  attributed  to  the  different  configuration  (1-exponential  2-compartment  vs. 
1-compartment  2-exponential)  and  data  set  (800  vs.  >  1700  man-dives)  used  to  estab¬ 
lish  the  MEPC  model  and  model  5  in  (6),  respectively. 

Finally,  the  method  of  maximum  likelihood  is  used  to  compare  the  predicted  risk 
of  DCS  among  different  decompression  tables.  Figure  6  shows  a  comparison  of  the 
probabilities  of  DCS  predicted  by  the  NLSC  model  for  standard  air  diving  to  150  fsw 
(4.54  ATA)  based  on  the  DCIEM  1983  air  decompression  tables  (12),  the  Kidd-Stubbs 
197 1  air  decompression  tables  (II),  and  the  U SN  standard  air  tables  ( 14)  as  a  function 
of  the  bottom  time.  (Table  2  fists  the  decompression  times  as  a  function  of  bottom 
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Depth  (fsvr) 


Fig.  5.  Probability  of  incidence  of  DCS  as  a  function  of  bottom  lime  for  standard  air  dives  to  150  fsw 
from  the  DCIEM  1983  tables  showing  the  difference  between  the  NLSC  model  ( solid  line),  the  MEPC 
model  ( broken  line),  and  model  5  in  (6)  ( dotted  line)  predictions. 

time  from  all  3  tables.)  All  decompression  procedures  produce  similar  predictions  for 
bottom  times  up  to  30  min.  The  Kidd-Stubbs  1971  decompression  times  are  more 
conservative  than  the  DCIEM  1983  tables  for  bottom  times  exceeding  45  min.  but 
less  conservative  in  the  moderate  bottom  time  range.  Therefore,  one  would  expect  a 
smaller  PD cs  for  the  longer  bottom  times  and  a  larger  Poes  for  the  moderate  bottom 
times.  This  predicted  behavior  is  shown  in  Fig.  6;  however,  the  differences  are  not 
very  large.  Similarly,  the  predicted  />DCS  using  the  DCIEM  1983  tables  are  lower  than 
that  using  the  USN  tables,  both  by  this  analysis  and  (hat  of  (15).  This  is  consistent 
with  the  longer,  and  more  conservative,  decompression  times  of  the  DCIEM  1983 
tables. 

A  powerful  application  of  the  method  of  maximum  likelihood  is  the  calculation  of 
decompression  tables  of  equal  PD cs-  Weathersby  et  al.  (6, 18)  have  used  this  method 
to  establish  a  set  of  air  decompression  tables  based  on  model  5  in  (6).  A  similar  set 
could  be  developed  for  any  model,  including  the  NLSC  model,  although  this  is  beyond 
the  scope  of  the  present  study. 

SUMMARY  AND  CONCLUSION 

The  method  of  maximum  likelihood  was  applied  to  compare  two  different  gas 
exchange  models  used  to  establish  safe  decompression  procedures.  These  models 
were  the  NLSC  and  MEPC  based  on  compartments  in  series  and  in  parallel,  respec¬ 
tively.  The  decompression  data  that  were  used  consisted  of  800  man-dives  with  21 
incidents  of  DCS  and  6  occurrences  of  marginal  symptoms  of  DCS.  A  comparison  of 


USE  OF  THE  MAXIMUM  LIKELIHOOD  METHOD 


311 


80 


ives  tuirO  fsw 
.vl.  thJMEPC 


:di/tions  for 
c^’are  fiore 

abut 
ct  a 
tom 

ices  a/e  not 

e  lowar  than 
s  coipistent 
)CIEM,  1983 


deflation  of 
this  method 
\  ^im|ar  set 


isie> 


liffej 

Fent 

gas 

hes^ 

m< 

[dels 

ill  el. 

re: 

Ipec- 

iV£K 

wi 

h  21 

>mpi 

\phon  of 

Bottom  Time  (min) 


Fig.  6.  Probability  of  incidence  of  DCS  predicted  by  the  NLSC  model  as  a  function  of  bottom  time  for 
standard  air  dives  to  150  fsw  from  the  DCIEM  1983  tables  (solid  line),  the  Kidd-Stubbs  1971  tables  [broken 
line),  and  the  USN  Standard  Air  Tables  ( dotted  line). 


TABLE  2 

Comparison  of  Total  Decompression  Times  for  Standard  Air  Dives 

at  150  FSW 


Bottom  Time,  min 

Total  Decompression  Time,  min 

DCIEM  1983 

KS-1971 

USN 

7 

2.5 

6 

3.5 

10 

9 

13 

3.5 

15 

18 

23 

5.5 

20 

25 

31 

11.5 

25 

43 

38 

23.5 

30 

57 

48 

34.5 

40 

88 

74 

59.5 

50 

128 

169 

88.5 

60 

178 

238 

112.5 

70 

228 

291 

146.5 

80 

271 

355 

173.5 

the  NLSC  and  MEPC  model  predictions  for  DCS  showed  no  significant  difference 
between  the  two  and  so  either  model  appears  equally  suitable  for  predicting  DCS 
with  this  limited  data  set.  Use  of  a  larger  data  set  would  improve  the  predictions. 
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allow  for  more  optimized  parameters,  and  perhaps  statistically  distinguish  the  pre¬ 
dictive  capability  of  different  models. 

A  striking  result  with  the  use  of  the  maximum  likelihood  method  is  the  difference 
in  the  predicted  probabilities  of  DCS  between  the  2  models  for  no-decompression  air 
dives  in  the  range  of  50-180  fsw  (1.51-5.45  ATA)  and  short-bottom-time  air  dives  to 
150  fsw  (4.54  ATA).  Because  of  the  nonlinear  gas  kinetics  of  the  NLSC  model,  a 
lower  probability  of  DCS  is  predicted  for  these  dives  with  this  model  than  with  the 
MEPC  model.  Subsequent  testing  indicated  consistency  with  actual  outcome  for  the 
NLSC  model  applied  to  the  no-decompression  air  dives  and  for  both  models  applied 
to  the  short-bottom-time  air  dives.  Further  studies  must  be  conducted  to  test  differ¬ 
ences  in  prediction  between  the  2  models  and  to  validate  the  extrapolation  of  model 
parameter  values  to  dives  that  fall  outside  the  range  of  the  data  set  used  to  establish 
their  values,  as  done  in  this  case. 

The  itithora  gratefully  acknowledge  the  advice  of  Dr.  B  C.  Eatock  in  the  preparation  of  this  paper.  DCIEM  Publication 
88-P-20. — Manuscript  received  February  1988;  accepted  May  1988. 


Tikusis  P,  Nishi  RY,  Weathersby  PK.  Emploi  de  la  methode  de  probability  maximale  dans 
I'analyse  de  piongees  a  1‘air  en  caisson.  Undersea  Biomed  Res  1988;  15(4):30 1  — 3 1 3. — La 
methode  de  probability  maximale  fut  utilisee  pour  t valuer  le  risque  de  maladie  de  decompres¬ 
sion  (MDC)  pour  un  choix  de  piongees  a  Pair  en  caisson.  Les  parametres  de  deux  modeles 
mathematiques  pour  prddire  la  MDC  furent  optimalis6s  jusqu’ft  I'obtention  du  meilleur  assen- 
timent  (tel  que  mesurd  par  la  probability  maximale)  correspondant  aux  incidents  de  MDC 
observes  dans  une  serie  de  piongees.  Les  donnees  de  decompression  employees  provenerenl 
de  800  piongees  humaines  contenant  21  incidents  de  MDC  et  6  cas  de  symptomes  marginaux. 
Le  premier  moddle  etudie  etait  base  sur  ie  principe  des  echanges  gazeux  mono-exponentiels 
dans  2  compartiments  montes  en  paralieie.  La  reussite  statistique  globale  k  declare  les  800 
piongees  humaines  fut  tres  similaire  pour  les  deux  modeles.  Les  predictions  de  security  pour 
les  piongees  ne  faisant  pas  partie  des  donnees  originates  etaient  differentes  de  celles  des 
mov'Mes  a  cause  des  variations  dans  la  cinetique  des  gaz.  Pour  les  piongees  courtes,  sans 
decon  oression,  les  compartiments  montes  en  serie  predisent  une  meins  grande  incidence  de 
MDC.  Ces  predictions  etaient  plus  consistantes  (s'accordaient  mieux)  avec  les  essais  subse- 
quents  que  les  predictions  avec  le  moddle  des  compartiments  en  parallelc.  Les  predictions  du 
moddle  avec  compartiment  en  serie  etaient  aussi  similaires  4  celles  obtenues  avec  le  moddle 
&  compartiment  simple  qui  fut  dvalud  par  la  Marine  amdricaine  avec  plus  de  1700  piongees 
humaines. 


Tikuisis  P,  Nishi  RY,  Weathersby  PK.  Empleo  del  metodo  de  probabilidad  maxima  en  el 
analisis  de  inmersiones  de  camara  con  aire.  Undersea  Biomed  Res  1988;  I5(4):301— 3 13. — Se 
empleo  el  metodo  de  probabilidad  maxima  para  evaluar  el  riesgo  de  enfermedad  por  des- 
compresion  (EPD)  en  inmersiones  de  camara  con  aire  electivas.  Se  optimizo  los  parametros 
de  dos  modelos  matematicos  para  predecir  la  EPD,  hasta  que  se  alcanzo  el  mejor  criterio 
(medido  por  la  probabilidad  maxima)  correspondiente  a  los  incidentes  por  EPD  observados 
en  una  serie  de  inmersiones.  Los  datos  de  descompresion  empleados  consistian  de  800 
inmersiones  de  humanos,  con  21  casos  de  EPD  y  la  incidencia  de  6  con  sintomas  marginales. 
El  primer  modelo  estudiado  se  basaba  en  un  intercambio  gaseoso  no  lineal  en  un  arreglo  en 
serie  de  cuatro  compartimientos.  El  segundo  modelo  se  basaba  en  un  intercambio  gaseoso 
monocxponcncial,  en  un  arreglo  en  paralclo  de  dos  compartimientos.  El  exito  estadistico 
global  para  describir  las  inmersiones  de  800  humanos  fue  similar  para  ambos  modelos.  Las 
predicciones  de  seguridad  para  inmersiones  que  no  formaban  parte  de  los  datos  originates, 
fue  distinta  para  los  dos  modelos  debido  a  la  diferencia  de  la  cinetica  de  los  gases.  La  serie 
de  arreglo  en  compartimientos  predijo  una  incidencia  menor  de  EPD.  para  inmersiones  de  no 
descompresibn  cortas.  Estas  predicciones  fueron  mas  consistentes  con  el  resultado  de  pruebas 
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subsecuentes,  que  las  del  modelo  de  compartimientos  en  paralelo.  Las  predicciones  del  modelo 
de  arreglo  en  serie,  tambien  resultaron  similares  a  las  del  modelo  de  un  compartimiento 
sencillo,  dos  exponential,  que  se  evaluo  con  mas  de  1700  inmersiones  de  humanos  por  la 
V.S.  Navy. 
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